Self-consistent calculation of the single particle scattering rate in high Tc cuprates 
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The linear temperature dependence of the resistivity above the optimal doping is a longstanding 
problem in the field of high temperature superconductivity in cuprates. In this paper, we investigate 
the effect of gauge fluctuations on the momentum relaxation time and the transport scattering rate 
within the slave boson method. We use a more general slave treatment to resolve the ambiguity 
of decomposing the Heisenberg exchange term. We conclude that this term should be decomposed 
only in the Cooper channel. This results in the spinon mass inversely proportional to the doping. 
It is showed that solving the equation for the transport scattering rate self-consistently, we find a 
crossover temperature above which we obtain the linear temperature dependence of the electrical 
resistivity as well as the single particle scattering rate. It is also shown that this linear temperature 
dependence of the scattering rate in the pseudogap region explains the existence of the Fermi arcs 
with a length that linearly varies with temperature. 
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It has been emphasized by many authors that 
the physics of high temperature superconductivity in 
cuprates [ij should be viewed as a doped Mott insula- 
tor 0-01 • This proximity to the Mott insulating phase 
at half filling, sheds some light on the underlying pairing 
mechanism of the superconductivity. One of the useful 
and rather successful treatments of the doped Mott insu- 
lator is the slave particle method [1-0|- It is well-known 
that this method results in the emergent internal gauge 
field that is strongly coupled to slave particles. Indeed 
it is the dressed slave particles by gauge bosons that is 
physical and observable. 

Since slave particles are strongly interacting with 
gauge bosons, the effect of gauge interactions on the 
physical properties of electrons should be taken into con- 
sideration in any realistic model. This effect has been 
extensively studied in the past three decades. For exam- 
ple, loffe and Larkin Q have studied how the physical 
properties of electrons, e.g. their electrical conductiv- 
ity, are related to the corresponding properties of slave 
particle, in the presence of the gauge field. Also, many 
authors have studied the effect of gauge fluctuations 
on the transport properties of electrons. Here we follow 
their method but with two modifications. First of all, we 
argue that the mass of spinous is inversely proportional 
to the doping, instead of being almost independent of it. 
As a result, the Fermi velocity of electrons should scale 
with doping outside the superconducting phase. The sec- 
ond difference is that we solve equations self-consistently. 
We show that either a simple scaling argument or exact 
numerical calculation, results in the transport scattering 
time that is linear in T above a crossover temperature, 
T* that scales with the doping. 
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From the Fermi liquid theory we expect dependence 
for the electrical resistivity. This reflects the stability 
of this phase at low temperatures. On the other hand, 
in the normal state above the optimal doping (strange 
metal phase of the hole doped cuprates), the resistivity 
exhibits a linear dependence on temperature instead of 
the expected parabolic dependence [ll. M [Toj . For temper- 
atures comparable to the Debye temperature and higher, 
this linear dependence can be explained through elec- 
tron phonon interaction. The Debye temperature is sev- 
eral hundred Kelvin while this behavior remains down 
to much smaller temperatures. So this is a fundamen- 
tal issue that cannot be explained only through simple 
electron phonon interaction. In slave boson approach, 
holons and spinons are strongly interacting with the in- 
ternal gauge field. One popular notion is that this inter- 
action can explain the desired behavior. By calculating 
the self energy of electrons due to scattering off the gauge 
field at finite temperature, we can compute the transport 
lifetime, as well as the momentum and energy relaxation 
time of quasiparticles resistivity. 

Some authors have studied the effect of gauge fluc- 
tuations on the transport properties of cuprates, and 
they find T^/^ temperature dependence for the resistivity 
[HQ. Here we closely follow their method, but with two 
important differences. First of all they assume the mass 
of spinons has a weak dependence on the doping and 
is inversely proportional to the exchange energy J. Us- 
ing the more general Anderson-Zou slave boson method 

0, we obtain a different behavior. In this method, 
the spinon mass is inversely proportional to the doping. 
This for instance means that the Fermi velocity outside 
the superconducting phase e.g. in the Fermi liquid phase, 
is proportional to the doping. The doping dependence of 
the spinon mass should be observed in the density of 
states at Fermi level. The second difference is that the 
problem should be solved self-consistently. In the rest 
of this paper, we have provided the details of calcula- 
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FIG. 1: Scattering rate. — Self-consistent calculation of the 
transport scattering rate 74^. and momentum relaxation scat- 
tering rate 7 versus temperature in the absence of external 
magnetic field. We assume that Cooper pairs are killed so 
we can study the crossover between the Fermi liquid and the 
strange metal phases. At low enough temperatures, holon 
gas condenses and we obtain FL phase. At high enough tem- 
peratures, we obtain linear temperature dependence for 'jtr 
and 7 signaling the strange metal phase, (a) 



verse transport lifetime 7*^ 
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ature is in units of Kelvin. Below the crossover temperature 
T* ~ 100-R', holons undergo Bose Einstein condensation and 
we obtain Fermi liquid behavior, 7tr ~ T^. Above T* , we 
obtain marginal non-Fermi liquid with 7tr ~ T. Interestingly 
T* scales with doping x, like the BEC transition temperature 
of the holon gas. 
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tion and we have presented an integral expression for the 
transport scattering rate ^tr that depends on itself. If 
we keep jtr inside the integral, along with the assump- 
tion 1/m* ^ X, we obtain a temperature scale T* above 
which we find linear temperature resistivity. Surprisingly 
T* scales with doping x, the same behavior that we ex- 
pect for the crossover temperature between the Fermi 
liquid and strange metal phases. Below this tempera- 
ture scale, holons condense and we recover the results of 
the Fermi liquid theory llj. We have solved the self- 



consistent equation for jtr and we verify behavior for 
the resistivity at low temperatures in the Fermi liquid 
phase and linear temperature resistivity above T* in the 
strange metal phase. Fig. 1 summarizes our result. 



III. METHOD 

Hubbard model is the simplest model that captures the 
physics of the Mott insulators. The Hamiltonian of the 
Hubbard model is defined as, 

The Hubbard model is defined as the following: 

H = U'^ Ui^^n^^l - t ^ c]^c^a + h.c. (1) 

Here (?, j) means site j is one of the nearest neighbors 
of site i. 

Now let us employ a more general slave boson approach 
to the Hubbard model [4|, 7]. In this approach an 
electron operator decomposes in the following way 



= A 



(2) 



where f^^ is spinon creation operator (assumed to be 

fermion), is holon creation operator and dj is dou- 
blon creation operator. We assume spinous are fermions 
and obey Fermi Dirac statistics and are electrically neu- 
tral, while holons and doublons are bosons and have +e 
and — e electric charges respectively. Spin a spinon cor- 
responds to a site with only one electron with spin up, 
holon represents an empty site and doublon corresponds 
to a doubly occupied site. As it is clear from the above 
definition, electron operator is invariant an under inter- 
nal compact U(l) gauge transformations {i.e. if we mul- 
tiply each creation operator by a local phase exp{i6i)), 
provided all slave particles carry the same internal gauge 
charge(i.e. 9i is the same for all slave particles at site i). 
Note that we must keep the bosonic statistics of holons 
and doublons after gauge transformation, otherwise we 
had SU(2) gauge freedom. In order to have fermionic 
statistics for electron operator, we have to implement 
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-d]di — 1 constraint at each site. 



On the other hand, doping implies {h\hi — d\di 
constraint. 

Since in the Mott insulator where onsite Coulomb re- 
pulsion is very large (U ^ t), the double occupancy 
is very costly, and these states can be removed in the 
low energy studies of cuprates. In other words, for large 
U /t limit of the Hubbard model, charge excitation gap 
is of order U and therefore we can integrate out dou- 
blons to obtain an effective action for spinon and holons. 
This process provides a systematic way to recover the t- J 
model and naturally reduces to the famous slave boson 
theory of the t-J model. 

In the t-J model, there is at most one electron at each 
site. To implement this constraint in the slave boson 
method, at we discussed before, empty states are treated 
as the charged bosonic particles, dubbed as holon. So we 
can take the non-doubly occupancy constraint by using 
the Lagrange multiplier method. The physical Hilbert 
space, contains three states at each site, occupied state 
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with spin up or down and unoccupied state (empty sites) . 
However, the Hilbert space of the slave boson method is 
much larger as we can have as many holons per site as 
we want and the constraint is implemented only in av- 
erage. Therefore the meanfield description of the slave 
boson method is incomplete and redundant. In the next 
section, it is shown that this redundancy is responsible 
for the emergent U{1) gauge field. Although at the be- 
ginning there is no kinetic term for the gauge field, upon 
renormalization and by integrating out the slave parti- 
cles, this term will be generated and the gauge field will 
have its own dynamics in that case. Slave particles are 
interacting strongly with this gauge field and their scat- 
tering off gauge potential gives them a finite lifetime. The 
single particle scattering rate is computed in section VI 
and we show that it is results in a transport scattering 
rate that is linear in temperature. 

One ambiguity in the t-J model is how to decompose 
the exchange term JSi.Sj, where Si is the spin operator. 
Most authors, decompose this term in direct and Cooper 
channels symmetrically. Within meanfield approxima- 
tion, they rewrite JJ2i j ^i-Sj as 



-VsjJ2{xQflM. + H.c. 

-3/8jJ2{AQafljl_^ + H.c. 
-3/8 J ^ 



(3) 



where x{j = T,a flafj,'^ and A{j = (^f-a,Jj^a- At 
small dopings, the above hopping implies 1/m* ^ Jxi-j 
which has a weak dependence on the doping. However, 
if we use the more general slave boson approach, it re- 
solves the above mentioned ambiguity and we have only 
one choice. In this formalism, the exchange term only 
decomposes in the Cooper channel. Using the Eq. [2], 
within the meanfield approximation, the Hubbard model 
can be rewritten as 

H = J2 udU, - 1 E (^^..4. + + H-^-) (4) 

in which we have used these notations Xij — h\hj — d\dj 
and A^^ = dihj + hidj. By integrating out doublons, 
it can be shown that (xb) ~ x, and (Afc) ~ t/U. 
A^jA^j represents the exchange term and after inte- 
grating out doublons it decouples only in the Cooper 
channel (we replace the exchange term by (A^ j) Af^- ^ 
jjP P + H.c. form). Therefore the hopping term for 
spinons is — < {xb) Yli j fi,afj,<^ and as results they have a 
mass 1/m* ^ 2tx. As we approach the half filling (un- 
doped material), the effective mass of spinons diverges 
signalling the metal-Mott insulator phase transition. 

So we end up at the the t-J model starting from the 
large U limit of the Hubbard model. It is believed 



that the t-J model captures the essential physics of the 
strongly correlated systems. This model is defined as the 
following 



t-J 



(5) 



where Pq is the Gutz wilier projection operator that 
removes doubly occupied states. Within the slave boson 
formalism, after removing doubly occupied states, elec- 
trons can be decomposed as c\ ^ — fj^hi along with the 

physical constraint on each site: hlhi + fj rrfi.a = 1 
which implements the Gutzwiller projection. Whenever 
cj ^ acts on an empty site, it annihilates one holon and 
creates a spinon with spin a. We cannot act further on 
the resulting state by cj since this operator has to kill 
a holon, but there is no holon anymore at that site. If 
we act Ci^a on a site that contains a spinon with spin cr, 
the operators annihilates the spinon and creates a holon 
at that site. So by acting projected electron operator we 
always annihilate one type of slave particle and create an- 
other one and therefore the number of slave particles at 
each site is conserved. Now we can rewrite the t-J model 
in terms of the new slave particles. Within meanfield 
approximation and by using Hubbard-Stratonovic trans- 
formation, we can decouple spinons (spin sector) from 
holons (charge sector) and we obtain the following effec- 
tive Hamiltonians for each sector 



Hh = - E tXsh\hj -^^iihh\hi (6) 

<i,j> i 
Hs = - E *Xhfl„fj,a -^fJ-sfl^Aa 

- E {J/2)A,i^,J)(fyl-flJl^+H.c.,{7) 



<l;J> 



where the following notations have been used 



Xh = [hlgK 



(8) 



(9) 



^^^^j)-l{fyi-fhfii) (10) 

At low temperatures, most of holons occupy the 
groundstate with momentum k = 0, therefore Xh ~ x. 
This model has been extensively studied in the literature 
and it is well known that this model leads to the d-wave 
pairing symmetry for spinons [3j|, i.e. As (±i) = As and 
As {±y) = -As. 



IV. SINGLE PARTICLE SCATTERING RATE 

Now, we discuss the effect of gauge fluctuations on 
the transport properties of cuprates in the strange metal 
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phase. Here we closely follow the method used by Lee 
and Nagaosa and a more recent approach by Senthil 
and Lee 

In the continuum approximation, the hopping part of 
the action of spinous and holons are written as 



dxdtfl {x) {^idt - VV2™: - /i.) [x) 
+ / dxdth^ {x) {-idt - V^/'^ml - fis) h (x) (11) 



As we pointed out, slave boson formalism has U(l) 
gauge theory. Under the gauge transformation, fi_„ -> 

e'^-/,,. and ^ e'^'h,. So xt:, ^ e'^'^-'^^^xtj ■ If we 

define ag{i) as Xij = \Xi.j\^^'^^''^\ where 6 — Ri ~ Rj, 
we have ag{i) — >■ ag{i) + 9j — 9i. This implies that in 
the continuum model, (x) — (a;) — 8^9 (x). There- 
fore we should add modify the continuum model in the 
following way to obtain a gauge invariant model 



J dxdtjl {x) {-iDt ~ Dj/2m: - fa {x) 
+ j dxdth'' (x) {-iDt - Df/2ml - h (x) (12) 

is the internal 



where = — 



ieintafj./c, where Cint 
gauge charge. We scale so that Cmt — e. 

To obtain an effective action for gauge particles we can 
integrate out spinous and holons. By expanding Z?^ in 
terms of and a^, it is clear that in the continuum 
model, the vector potential is minimally coupled to mat- 
ter field, i.e. the gauge field is coupled to the current 
carried by quasiparticles and we have the following inter- 
action between gauge field and spinous 



dxdtJ!^ (x, t) {x, t) 



(13) 



where = (0, ax, ay) and = (n, J^, Jy), where = 
P f, Js — ip (v/2ms^ f + H.c. . We have a similar 
term for holons. Now we use the gauge freedom to set 
= and choose the Coulomb gauge V.a = 0. Up to the 
second order perturbation theory the above interaction 
terms leads to an following effective action for the gauge 
field proportional to the following form, 



dxdrdx'dr'ai {x, r) ( J* (x, r) , J-' {x' , r')) aj {x, t)(14) 



This is only the contribution coming from paramag- 
netic response. There is a similar contribution from 
the diamagnetic part. Lee and Nagaosa have stud- 
ied this problem in detail. Following their method, 
we can find the gauge field propagator -D^^^ (r, r) = 
(T^a^ (r, r) (?-', 0)). it can be shown that the 
gauge field propagator can be written as Dij (q) = 
(n^ + n,i)^r.\ where 

ir^ = - {Tr [Jl (r, r) (0, 0) - S^.n^S (r) S (r)] ) (15) 



and similarly for H]^ . In the Coulomb gauge the spatial 
part of the gauge field is transverse and can be written 
as 

{q,i^) = - qiqj/q^) d'^ {q,uj) (16) 

where we have defined {q, ut) as 

D^{q,u;)^[uUq,u;) + Ul{q,u;)]-' (17) 

For small q and w we can use the following approximation 

nf (g, u) = iLua^ [q, uj) - x?q^ - Ps^c (T) (18) 

where aJ is the transverse conductivity and x? is the 
Landau diamagnetic susceptibility of the spinon system, 
which equals l/247rms in 2D fcrmion systems, ps^c is the 
condensation fraction of the spinon system in the paired 
state and in the normal state it is zero. Similarly for 
holons we have 

nf {q, oj) = i H (Jh (?) ^) - Xhq^ - mhPh.c (t) (19) 

where Xh ~ ^ (0) /487rms in which n (0) is the occu- 
pation number of the ground state and ph^c is the su- 
perfiuid density of holons, which is zero in the strange 
metal phase as well as the pseudogap phase. We see that 
the propagator of the gauge field that determines the 
transport scattering rate of spinous, depends on the op- 
tical conductivity of spinous. On the other hand, optical 
conductivity itself depends on the transport scattering 



rate, since we have as{q,uj) — v^kp/ y 'ij^^ + v'^q'^. If 

we neglect 7tr in the as {q,uj), we recover the results of 
the Ref. [1], except that since the spinon mass is much 
larger in our case, the crossover temperature is smaller. 
In Ref. 9] , Lee and Nagaosa have calculated the scatter- 
ing rate of holons due to interaction with the gauge field 
and they report a linear temperature behavior above the 
BEC transition. They have also studied the transport 
scattering rate of spinous and they report T*/"^ in the 
strange metal phase, which is slightly different from the 
linear temperature behavior. Now let us calculate the 
self energy of spinous due to their interaction with the 
gauge field. Then by looking at its imaginary part we 
can deduce the scattering rate. According to Senthil and 
Lee, up to the first order approximation, the self energy 
of spinous due to their scattering off the gauge field is 
described by the following expression 



T,"{K,n = 0) = -Try dujdq{v^ x qy D" {q,uj) 
xA{K-q, -uj) {n,, {uj) + n,^ (^)) (20) 

where n^^ and n^j^ are Bose Einstein and Fermi-Dirac 
distribution functions respectively. Near the Fermi sur- 
face the above integral reduces to 



7 = TTVpK^ 



dw- 



1 



sinh (3uj 



dqD" {q,uj) 



(21) 
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where D" {q,uj 
a = as + a/j 



Xh and Pc = Pc, 
compute we should multiply D" {q,u}) by ■ So 



in which 
Pc,h- To 

2 



we have 



Itr 



duj- 



1 



sinh jiuj 



dq 



q 



D" (g,^)(22) 



Let us assume that ^ Kp/m* ~ 2txKp and at 
small dopings, ~ tt. When g <C ^tr/Vp, CTs = 
Vpkp/^tr- On the other hand, in the absence of holon 
and spinon condensation, wcr(g,aj) ^ X^'f': so q ^ 

/ \l/2 

i^i? ( 2^ j • Therefore the q < q* ^ Itr/Vp region con- 

1 /2 

tribution to the integral over q is of order {"ftr/Vp ) 
For g* comparable to a fraction of tt this is the leading or- 
der contribution. This corresponds to 7*'' (T) cx oc te. 
As long as this is satisfied we have 



3/2 

7tr Of 'r^'f 



da;- 



;l/2 



sinh (3(jj 



(23) 



A simple scaling argument verifies that ^tr ^ T sat- 
isfies the above self-consistent equation for 74,.. The 
crossover temperature can be found from 7*'' (T*) ^ 
Vp oc tx. Therefore T* cx tx, which is the expected be- 
havior for the crossover temperature between the Fermi 
liquid phase and the strange metal phase. 

Numerically, we solved this self-consistent equation ex- 
actly and we again we obtained linear T dependence of 
7ir above a temperature scale T* [x) that scales with 
doping X (see Fig. 1). On the other hand the BEG 
transition temperature of holons scales with doping as 
well. In summary, above a temperature scale compara- 
ble to the Bose Einstein transition temperature of holon 
gas, we obtain linear temperature transport scattering 
rate. It should be mentioned that the scattering rate of 
holons is linear in T even in the previous calculations 9] . 
To obtain the physical quantities from the corresponding 
quantities for spinous and holons, we should recombine 
them in a particular way, that is called loffe-Larkin for- 
mula. This formula tells us that the physical conductivity 
is related to that of holons and spinous in the following 
way (T^^ = (77^ + '^h^- conductivity we have 

a = ne^Ttr/ra and Tt,- = 1/jtr- Since the scattering rate 
of holons and spinous are both linear in T, the dc con- 
ductivity is is linear in temperature as well. In Fig. 1 we 
have presented our numerical results. 

On the other hand, below the BEG transition tem- 
perature, Pc,/t 7^ 0- Numerically, we obtain behavior 
for the jtr in this phase, which is the right sign and is 
expected for the Fermi liquid phase. 
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FIG. 2: Length of the Fermi arc. — Length of tlie Fermi 
arc. In the pseudogap region, we obtain disjoint Fermi seg- 
ments, at which excitations are gapless. The length of these 
segments increases with temperature and the scattering rate 
7. 1(a) [ The length of the Fermi arc vs. normalized tempera- 
ture (T/T*). The blue line our theoretical calculation and red 
triangles are experimental data taken from Ref. Below 
the superconducting transition temperature Tc, the scattering 
rate 7 is very small due to the condensation of quasiparticles. 
Therefore only very near nodal points we have gapless ex- 
citation. Above Tc, there is macroscopic condensation and 
the scattering rate varies linearly with temperature. As 7 
increases, the length of the Fermi arc increases. At T* , we 
obtain closed Fermi surface, (b) The length of the Fermi arc 



vs. scattering rate 7. Below the superconducting transition 
temperature Tc, the scattering rate 7 is very small due to 
the condensation of quasiparticles. Therefore only very near 
nodal points we have gapless excitation. Above Tc, there is 
macroscopic condensation and the scattering rate varies lin- 
early with temperature. As 7 increases, the length of the 
Fermi arc increases. At T* , we obtain closed Fermi surface. 



V. OBSERVATION OF FERMI ARCS 

In the underdoped cuprates and above the supercon- 
ducting region, the state of matter is very unusual. In 
one hand, it is a metallic phase, confirmed by transport 
experiments. On the other hand, there is no closed Fermi 
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surface. Along some directions and segments of the Fermi 
surface, excitations are gap less, while on other directions, 
excitations are gapped |12l - ll4| . c-axis transport measure- 
ments also show a gap in the excitations that implies 
a bound state of Fermions. It has been discussed by 
Norman et al [13] that if the scattering rate of quasi- 
particles 7 varies linearly with temperature, one can ex- 
plain the observations. As we discussed in the previous 
paragraph, gauge fluctuations can cause such a behav- 
ior. Some authors [14] have discussed that the electron 
phonon interaction may play an important role here and 
result in the linear T dependence of 7 above some at high 
enough temperatures. Here we show that the interaction 
of slave particles with gauge bosons and the scattering 
of the quasiparticles from the d-wave potential results in 
such an exotic behavior. Using the method introduced 
in the previous section, we can compute the scattering 
rate of electrons in the superconducting region as well as 
the pseudogap phase. In the superconducting phase, the 
scattering rate is very small due to the condensation of 
both spinous and holons. In the pseudogap phase how- 
ever, it is comparable to the pseudogap, varies linearly 
with T and cannot be neglected. In this region we assume 
that there is a local pairing potential that electrons scat- 
ter off. This assumption leads to the following expression 
for the electrons Green's function: 



G{k,u:)-^ =oj-ek + ij , \ . (24) 

where, 7 (T) is the scattering rate at temperature T, en- 
ergy of free quasiparticles. We have assumed d-wave pair- 
ing, i.e. Afe (r) = A (T) (cos (k^) — cos {ky)). Using the 
above expression we can compute the spectral function 
and from the position of its peak we can read the en- 
ergy of interacting quasiparticles. When = 0, it can 



be shown that as long as 7 > -s/SAfc (T), the maximum 
is at = and as a result we have gapless excitations 
along that direction. At nodal points A^ — and we 
can always satisfy this equation at that point. Since 
7 (T) ^ A (T) in the superconducting phase, this condi- 
tion is only satisfied very near nodal points and the length 
of Fermi arcs is negligible. If we parameterize A^ by 
the angle, we have A (0, T) = A (T) cos (20), and nodal 
points are located a.t (p — {±7r/4, ±37r/4}. Therefore the 
length of arc is L (T) = 8sin~^ (7 (r)/%/3A(T)). The 
crossover temperature of the pseudogap phase, T*, can 
be found by solving 7 (T*) = ^/3A (T*). Fig. 2 summa- 
rizes our results. 

In conclusion, we reexamined the effect of gauge fluc- 
tuations on the single particle scattering rate. Using the 
mass of spinous proportional to the doping and doing a 
self-consistent calculation for the self-energy of quasipar- 
ticles scattering off the gauge field, we found the linear 
temperature dependent transport scattering rate of elec- 
trons above a crossover temperature that scales with the 
doping. Below this crossover temperature, we obtained 
behavior for the transport scattering rate that recov- 
ers the Fermi liquid behavior. We also found a scattering 
rate that is linear in T above the crossover temperature. 
We showed that the linear dependence of the scattering 
rate explains the existence of the Fermi arcs in the pseu- 
dogap phase. The success of our model emphasizes on 
the importance of including gauge fluctuations in under- 
standing the underlying physics of cuprates. We predict 
that the Fermi velocity in the Fermi liquid phase varies 
linearly with the doping. 
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